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We examine timing noise in both magnetars and regular pulsars, and find that there exists a component of 
the timing noise ((Txn) with strong magnetic field dependence (cttn ~ BlSlT^'-) above B„ ~ 10'--^G. The 
dependence of the timing noise floor on the magnetic field is also reflected in the smallest observable glitch size. 
We find that magnetospheric torque variation cannot explain this component of timing noise. We calculate the 
moment of inertia of the magnetic field outside of a neutron star and show that this timing noise component may 
be due to variation of this moment of inertia, and could be evidence of rapid global magnetospheric variability. 



Introduction. — The spin evolution of pulsars and magne- 
tars is studied by long term monitoring of pulse arrival times, 
allowing the spin period, spin-down rate, orbital motion and 
astrometric variations of the neutron star to be inferred. The 
variation of these pulse arrival times from the best fit models 
of the timing evolution is referred to as timing noise. 

Timing noise was first observed in the Crab pulsar 1 1 1 
and has been identified as a ubiquitous property of pulsars 
l2|. Recently much effort has been made to understand the 
stochastic timing noise in millisecond pulsars, as timing ar- 
rays of low-noise recycled pulsars will allow the detection 
of nanohertz gravitational waves [3|. Millisecond pulsars 
typically have relatively low inferred surface magnetic field 
(Bo ~ lO' - lO^G), due to age and accretion history. 

Most isolated pulsars have spin-down inferred surface mag- 
netic fields in the range B,, ~ 10" - lO'^G, while the high- 
est known magnetic fields are possessed by magnetars ~ 
10'^* - lO'^G). Magnetars include Anomalous X-ray Pulsars 
(AXPs) whose X-ray luminosities dwarf their spin down lu- 
minosities and most likely originate from internal magnetic 
field decay |4|; and Soft Gamma Repeaters (SGRs), hard X- 
ray transient sources that can undergo extreme outbursts. 

Mode-changing and nulling behaviors in pulsars have been 
linked to spin-down torque variations f5^ '51, implying that 
such behavior is indicative of variations of the open field line 
regions of the magnetosphere. Recent simultaneous X-ray 
and radio observations of PSR B0943-I-10 [7 | have challenged 
existing pulsar emission theories and provided evidence that 
such variability is indicative of rapid global variations of the 
magnetosphere. 

In this Letter we study the timing noise for radio pulsars and 
magnetars, inferring a source of timing noise in high B-field 
pulsars and AXPs which coiTelates strongly with magnetic 
field. We examine physical models to explain the dependence 
of this timing noise on the spin-down inferred magnetic field, 
and identify this timing noise as evidence for global magneto- 
spheric variability. 

Timing Noise Analysis. — Long term pulsar timing irregu- 
larities have enjoyed a long history of detailed analysis [see 
e.g. [8 UT2II . Here we adopt the approach of modelling tim- 
ing noise as a random walk process lUl [13] [141 . In particular, 
we utilize the formalism of Cordes lfT4l in order to define the 
random walk strengths of various processes: 5pN = 



for a random walk in phase; 5fn = RyiS^)/, for a random 

walk in frequency; and Ssn - R ((^fi)^), for a random walk in 
spin-down rate, where R is the occurrence rate of the random 
walk steps, and { ■ ) indicates an ensemble average. The quan- 
tities 64>, 6Q., and 6Q. denote stochastic variations in phase, 
frequency, and spin-down rate, respectively. 

These strengths can be estimated from timing parameters 
of a given pulsar and represent the strength of the assumed 
random walk 

SpN - ^ClmCrl^T (1) 

5fn - nciy^^T-\ (2) 

5sN - l20Cl„a-',^T-\ (3) 

where, T is the time span over which the observations were 
(uniformly) taken, ctxn is the rms phase residual of the data for 
a given timing solution, and Co,,,,, Ci.,„ and C2,m are correction 
factors [14, 15 1 to compensate for the power of the random 
walk that has been removed by the mth order polynomial fit 
that occurs when the residuals are determined. 

Previous analyses [8, .10, 12] have shown that simple ran- 
dom walk processes cannot explain the totality of timing 
noise, with the different strengths sometimes varying strongly 
even for the same pulsar as a function of timespan T. We 
argue that if some random walk timing noise component be- 
comes dominant as the magnetic field increases 8 orders of 
magnitude from ~ lO^G to ~ lO'^G this should result in a 
lower bound in the distribution of the random walk strength 
that scales with the magnetic field. 

In Figure [T] we show the random walk strengths for vari- 
ous pulsars versus the spin-down infeiTed surface dipole mag- 
netic field strength. The strength of the frequency noise 
and spin-down noise has been normalized by and Cf, 
respectively, so that the strengths can be compared across 
pulsars with differing timing profiles. Here, we utilize pub- 
lished timing data from Jodrell Bank Observatory [8|, and the 
Parkes 64m radio telescope flBl - BTl . We also include the X- 
ray timing results from the well-timed AXPs: IE 1841-045 
11251 . RXS J170849.0-400910 [22|, 4U 0142-h6I|231, and IE 
2259.1+586 |24|. We do not include the timing for AXP IE 
1547.0-5408 |25l, as the timing observations have only been 
taken post-outburst, nor AXP IE 1048.1-5937 as the 
timing solution presented in the literature was fit using fifth or- 
der splines rather than a simple polynomial fit, due to instabil- 
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FIG. 1. Effective random walk strength for various pulsars as a func- 
tion of spin-down inferred magnetic field strength (Bo). The fre- 
quency noise (FN) and spin-down noise (SN) strengths have been 
normalized by Qr and Cf, respectively, to be comparable between 
systems. The blue points denote normal pulsars, with radio tim- 
ing information taken from Hobbs et al. (§], Yu et al. 1161 , and the 
Parkes multi-beam survey 1 17^- 211. The red circles label AXP timing 
epochs, with X-ray timing information taken from Dib et al. 1 22 , 23 1, 
and Gavriil and Kaspi |24|. The effective strength of the frequency 
noise random walk shows a strong trend with magnetic field strength. 
In particular the lower bound of the frequency noise strength distri- 
bution appears to rise sharply for Bg > 10'^^ G as (Sfn/^^^)™™ ~ S*, 

such that (0-TN)m,n ~ BIQ. 



ity of the spin-down of IE 1048.1-5937. We also ignore the 
timing properties of SGRs (unlike Shannon and Cordes 0), 
which are only timed in a phase connected fashion for short 
periods following outbursts where timing noise would likely 
be elevated by the burst activity. In quiescence, where the tim- 
ing noise would not be dominated by such events, SGRs are 
faint and have not been observed with sufficient regularity for 
a phase connected solution to emerge. 

While there is no obvious correlation between magnetic 
field strength and the random walk strengths for phase noise 
(5pn) or spin-down noise (Ssj^/Cf-) in Figure[T| there appears 
to be a weak correlation of the FN strength (Sp^^/Q.-^) with Bg 



across the range of magnetic fieldM We note, however, that 
above a field strength of Bg 2" 10 ^ G, the lower bound of 
the FN strength distribution rises sharply with (5FN/i^^)min ~ 
q2/q6 ~ B^. We focus on FN as the change in 5fn along this 
lower envelope is larger than the intrinsic scatter of the distri- 
bution for FN (particularly due to the inclusion of AXP timing 
data). Thus, we infer that timing noise has a component which 
depends strongly on the magnetic field and becomes dominant 
as the magnetic field increases such that (cr jN)min ~ bID.T^^^ 
above B^ > 10'^^ G. These scalings are within the 2cr con- 
fidence intervals for the timing noise scalings inferred by j51 
for magnetars, except for the dependence on T, which may be 
different due to their inclusion of SGRs. 

Caution must be used in interpreting such scalings as B^ ~ 
Q/Q^, and is only inferred from the spin observables. With 
this in mind, we consider two different physical models of 
timing noise due to magnetospheric variability, torque varia- 
tion, and moment of inertia variation. 

Magnetospheric Torque Variation. — Spin-down torque 
variation and mode-changing are associated with perturba- 
tions of the open field lines ||5]|6l. The effect of these torque 
variations on the timing noise has been explored in detail by 
Cheng 1 27 1, who showed that the frequency noise strength for 
torque variability should scale as ^PN/fi^ ~ O^/fi"' ~ B^Q.^, 
which is inconsistent with the scaling discussed above. Thus, 
while this source of noise may dominate in some pulsars, par- 
ticularly those where significant pulse-shape changes are ob- 
served, it is not the source of the high magnetic field timing 
noise floor evident in Figure [T] 

Magnetospheric Moment of Inertia. — The angular momen- 
tum content of a small volume of the magnetosphere in the 
inertial frame [28 1 is liL = (r xS/c^)dV where S is the Poynt- 
ing vector and dV is the volume element. In a co-rotating 
ideal-MHD magnetosphere with angular velocity Q, we can 
evaluate this as 



r X (£2 X r) 



4nc^ 



- (r X B) 



B ■ (Q X r) 

47rc2 



dV. (4) 



In cylindrical coordinates, (■U7,z,(p), the component of the 
angular momentum along the rotation axis is dL, = 
l^n7^Q(B|, -I- Bl)dV. Thus the angular momentum of the 
magnetosphere in the direction of the rotation is given by the 
integrated angular momentum of the magnetic energy density 
with the component of the magnetic field in the azimuthal di- 
rection removed, as the component of the magnetic field in the 
direction of motion should not contribute to the Poynting flux 
in that direction. 



' The weak correlation of this timing noise with magnetic field for lower 
field strengths has been noted previously, see e.g. Figure 9 of Hobbs et al. 
| 8 |, where timing noise is measured by the parameter o"j(10yr). The onset 
of a sharp increase in o"-(10yr) above ~ lO'^-^G can also be seen in this 
figure, but is not mentioned by the authors. In Figure^the trend is more 
clear due to the inclusion of the AXPs. 
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For a dipole with magnetic moment fig, where \fio\ - Bgr^^ 
and 2B„ is the field strength at the magnetic pole at the NS sur- 
face tns, we calculate the magnetospheric angular momentum 
(for light-cyUnder radius Wic = c/Q » rjsjs) using (Hh 
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'"NS 



-Qiln-[fXo-n]tio). 



(5) 



For the aligned rotator this gives LB.aiiened 



(16/15)Z?^r^j£l/c^, while for the oblique rotator we have 

LB,obMque - (17/15)B2r5^n/c2. 

For a magnetic field strength Bi5=Bo/10'^G, neutron star 
radius re = rNs/(10^cm) and mass Mi .4 = Mns/(1-4M0) we 
can compare the moment of inertia in the magnetosphere to 
that in the neutron star /b//ns - Ls/L^^s, 



'B, aligned 



(6) 



for the aligned rotator, where 77 = /ns/(A^ns'"ns) - 2/5 for a 
uniform rotating sphere. While this ratio depends strongly on 
the NS radius, typical equations of state which allow masses 
as large as the observed 2Mq, have radii varying by at most 
~ 20% over the range of expected NS masses 1291 l30l . 

Magnetospheric plasma has at least the Goldreich-Julian 
density pa = ■ B/[2nc(l - (m/mi^cf] EH. Near the 
NS surface, vr <K vj^c, the ratio of the Goldreich-Julian 
plasma energy density to the magnetic field energy density 
is £Gj/£mag ~ lO-^'^(n/s-^)(juJlO^^^Gcm^)-Um/l0^cm)\ 
The radius cTeq at which the energy density of the plasma is 
comparable with the magnetic field energy density can be es- 
timated as tzTeq = wlc Vl - ^■c'^mg / QioOP-e). For a pulsar with 
Bo ^ lO'^G, this gives tiTeq - 0.999tjTz,c which is extremely 
close to the light cylinder Relativistic force-free solutions of 
the magnetosphere that take into account the poloidal currents 
1321 lead to similar results as their differences are concentrated 
near the light cylinder, while most of the mass and energy 
density is near the NS surface. Thus the contribution of the 
Goldreich-Julian plasma to the moment of inertia of the mag- 
netosphere can be safely ignored. 

Magnetospheric Variability. — Mode-changing and nulling 
events are related to rapid variability of the open field line re- 
gion of the magnetosphere lHJISl. Recent observations have 
shown that such rapid variability may be a global magneto- 
spheric phenomenon |7|, and not simply confined to the open 
field line region. Thus, the magnetospheric moment of inertia 
could also vary on a short timescale. 

The equation of motion for rotation is ^ [/(f)Q(f)] - N(t), 
where A^(f) is the external torque on the neutron star. We de- 
fine I - Ic + 61{t) where SI{t) is a stochastic component to the 
moment of inertia that we will associate with magnetospheric 
variations, and Ic is the moment of inertia of the NS that is 
strongly coupled to the crust such that the coupling timescale 
is much shorter than the timescale of the variability. This 
strongly coupled component is the part of the star that can 
then respond eff'ectively to the moment of inertia variation. 
We also define Q.(t) - Q.s{t)+SQ.{t) where flj (f) is the smooth 
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FIG. 2. A measure of the moment of inertia variability A^'y^A/mis/Zc 
versus the surface magnetic field strength B„ for the same pulsars 
and magnetars shown in Figure [T] A'ly^ is the number of moment 
of inertia random walk steps that occur per year. We compare this 
to the fractional moment of inertia contained in the magnetosphere, 
h/Im (assuming an aligned dipole for a IAMq mass neutron star 
with radius 12 km). 



polynomial behavior of the angular velocity and the stochas- 
tic component is 6Q.(t). Ignoring spin-down torque variations 
that could occur from magnetospheric variation near the open 
field lines, we solve for the 6Q.(t) to first order in 61 /Ic and 
6Q./D.S, assuming that the torque takes the form = aQ" 
where n is the braking index of the neutron star. This gives 



5Q(f) = -Q5(0 



6I(t) 



j nnsit')'-^dt' 



(7) 



The second term is much smaller than the first term as long 
as the observed time span is short compared to T <«c = 
Qs /2Q.S > 10"* years for typical magnetars. 

Modeling the stochastic variation of the moment of inertia 
as a random walk we have 61{t) - Yjj AIjH(t - tj), where A/y 
and tj are random amplitudes and times, while H{t) is the unit 
step function. The first term of equation (|7]) then corresponds 
to a random walk in frequency, while the second term is a 
random walk in frequency derivative. Considering realistic 
observing spans the first term must dominate. Using the defi- 
nition of 5 FN, the variability of the moment of inertia can then 
be expressed in terms of the frequency noise random walk 
strength. 



yl(S^/D?)xlyr. 



(8) 



where we have chosen to use the parameter A^iy, = R x lyr, 
the number of random walk steps per year. This value is 
plotted against magnetic field strength in Figure |2] as well 
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as various fractions of the magnetospheric moment of inertia 
Ib/In5, assuming an aligned magnetic dipole and a neutron 
star mass I.4M0 and radius 12 km. We note here that above 
Bg ^ lO'^'^G the moment of inertia variability is bounded by 

■ 1 /2 

A^j^j. A/,-n,s/4 ~ (0.1 - 1)/b//ns- If the timing noise observed at 
this lower bound is due to variability of the moment of inertia 
it follows then that Nl!^^(MrmJh)(lNs/Ic) -0.1-1. 

Assuming A^iyr ~ 10^ - 10^ to reflect the variability 
timescales observed in pulsar nulling or mode-changing 15J, 
and the strongly coupled fraction of the NS moment of inertia 
to be /c//ns ~ 0.1|"27l over the variability timescale, we can 
estimate the rms amplitude of the magnetospheric moment of 
inertia variability, AI^^s ~ (lO""* - 10"'') /b. 

The magnetosphere is expected to be dynamic near the light 
cylinder, due to reconnection and instability [32, .331. This 
may also contribute to the timing noise through magneto- 
spheric moment of inertia variation, however the amplitude of 
this contribution to the random walk strength is suppressed a 
factor (rNsO/c)^ ~ lO^^s^ Q?, due to the much weaker contri- 
bution to the moment of inertia near the light cylinder. We find 
that this component is too small to contribute substantially to 
the timing noise discussed above. This implies that the vari- 
ability of the moment of inertia must instead occur near the 
NS surface, where the magnetospheric moment of inertia is 
largest. 

Comparison to Glitch Sizes. — Glitches in pulsars and mag- 
netars are impulsive increases in the rotation frequency. They 
are thought to represent a sudden transfer of angular momen- 
tum from a more quickly rotating super-fluid component to 
the NS crust. While glitches can vary across several orders of 
magnitude in size for a single pulsar, the smallest detectable 
glitch sizes can also serve as a measure of timing noise. 

The smallest detectable glitch size can be estimated from 
the timing noise level, by comparing the (pre-fit) phase change 
due to the stochastic noise over the data span required to in- 
ferred the existence of small glitch, to the phase change due to 
the ghtch itself, A^m(At) ^ S\!^{At)^l- 1 ^fU. < AOgHtchAf 
for frequency noise. If the expected phase change due to 
the noise is larger than the phase change due to the glitch, 
then a glitch cannot be definitively identified. Oversimplify- 
ing the details of observational cadence and analysis we can 
then estimate the smallest observable glitch size by assuming 
that several (~ 3) time-of-anival observations, with cadence 
~ 1 month, on either side of a small glitch are needed to 
characterize it. Combining these assumptions with the tim- 
ing noise levels estimated above we arrive at the estimate 
(AQ/Q)gii„h > (0.02 - 0.2)/b//ns. 

In Figure |3] we show the relative glitch sizes (AQ/Q) as a 
function of magnetic field for the glitching pulsars listed in 
the literature lfT6l l34ll . We also include glitches from AXPs 
4U 0142+61 fSTl, IE 2259.1+586 |24|, IE 1841-045 Il22l . 
RXS J 170849.0-4009 10 122J and IE 1048.1-5937 l26l. Also 
plotted is the ratio of the magnetospheric moment of inertia to 
the NS moment of inertia (/b//ns) as a function of B,,, as well 
as 10%, and 1% of this value. We find that for B,, > 10'^ G, the 
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FIG. 3. Relative glitch size, AQ/fi, versus the spin-down inferred 
magnetic field strength B„, for radio pulsar glitches, taken from Es- 
pinoza et al. OH and Yu e/ aZ. 1 1 6 1 , and for AXP glitches 1 22 , 24 26 1 . 
The ratio of the magnetospheric moment of inertia to the total mo- 
ment of inertia (/b//ns)> is also plotted as a function of B„, along with 
10% and 1% of Ib/Ins- The minimum observed glitch size serves as 
a proxy for the timing noise, scaling with (ASl/Si)^i„ ~ 0.3/b//ns for 
So > lO'^G and is consistent with our estimates of timing noise. 



minimum observed glitch is roughly given by (AQ/Q)giitch ^ 
0.3/b //ns, which is consistent with our estimates above. Thus, 
we find consistent evidence from glitches for an increase with 
Bg of the timing noise floor. 

While in principle a glitch due to change in the magneto- 
spheric moment of inertia could be detected, it would require 
a large (> 10%) change in the total magnetospheric moment of 
inertia to be above the minimum detectable glitch size thresh- 
old set by the timing noise. Such an event would almost cer- 
tainly be accompanied by torque variations and particle out- 
flows, as seen during giant flares in SGRs, which would dom- 
inate the timing change due to the magnetospheric moment of 
inertia. 

Discussion. — We have examined pulsar and AXP timing 
noise measurements reported in the literature and shown that 
a component of timing noise exists which depends on the spin- 
inferred dipole surface magnetic field strength (S m/Q.^ ~ B^ 
such that cTxN ~ BIQT^^^). This timing noise component be- 
gins to dominate at B„ ~ 10^--^ G, and is responsible for a 
sharp rise in the floor of the timing noise values across both 
pulsar and AXP populations. 

This provides yet another connection between high-B radio 
pulsars and AXPs, demonstrating a continuum of behaviors in 
these neutron stars, as independently suggested by, for exam- 
ple, quiescent X-ray luminosities of these objects ll35l . further 
'unifying' radio pulsars and AXPs l36\. 

Variations near the open field lines can lead to mode chang- 
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ing and torque variability f5\, however this would result in a 
frequency noise strength that scales as Spt^/Q.^ ~ B^Q?, and 
cannot explain the observed timing noise dependence. 

We have shown that the magnetospheric moment of inertia 
is /b ~ IO^^B^j/ns, and have proposed a model of magneto- 
spheric moment of inertia variation that fits the observations 
for both the observed timing noise strengths and the size of the 
smallest observable glitches in high magnetic field systems. 
This variation must occur near the neutron star surface, where 
the moment of inertia contribution is largest, as we find that 
timing noise due to variability confined near the light cylinder 
is too small to contribute significantly. Internal field evolution 
could also act as a source of moment of inertia variability. 

This timing noise dependence provides evidence of rapid 
variability throughout the magnetosphere through moment of 
inertia variations. By assuming a rate similar to the known 
variability in the open field line regions of some pulsars ||5l|6l 
we can estimate an amplitude for the variability of the magne- 
tospheric moment of inertia Al^ms ~ (10""^ - 10'')/b. 

Recent observations [71 have also provided evidence that 
there exists rapid global variability in pulsar magnetospheres. 
We suggest that rapid global magnetospheric variability, per- 
haps due to reconnection or variable currents which are known 
to exist from pulsar simulations ||33]| . acts as a source of tim- 
ing noise through moment of inertia variations. 
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